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Abstract 

A  new  mechanism  for  the  thermal  desorption  of  molecular  hydrogen 
from  the  monohydride  phase  on  Sl(100)  has  been  identified.  The  unusual 
first-order  desorption  kinetics  that  are  observed  are  due  to  the  irre¬ 
versible  excitation  of  a  hydrogen  adatom  into  a  delocalized,  two- 
dimensional  band  state  on  the  surface  with  an  activation  energy  of  47 
kcal/mol.  The  desorption  reaction  occurs  between  this  excited  hydrogen 
adatom  and  a  second,  localized  hydrogen  adatom.  The  mechanism  was 
verified  by  adsorption  of  atomic  hydrogen  into  the  band  state  and  the 
observation  of  reaction  with  localized  deuterium  adatoms  at  a  tem¬ 
perature  below  that  at  which  desorption  from  the  surface  normally 
occurs. 

PACS  numbers:  68.45  Da,  82.20  Pm,  82.65  My;  82.65  Yh. 


In  this  Letter,  we  examine  the  Isothermal  recomblnatlve  desorption 
of  hydrogen  from  SI (100)  and  Identify  a  new  mechanism  for  the  desorption 
reaction.  On  metal  surfaces,  recomblnatlve  desorption  proceeds  via 
migration  of  adatoms  from  site  to  site  until  two  adatoms  in  adjacent 
sites  ultimately  recombine  and  desorb.  We  shall  show  that  this  mecha¬ 
nism  is  inconsistent  with  the  isothermal  desorption  kinetics  of  hydrogen 
from  Si(100),  measured  between  685  and  740  K.  We  propose  and  verify 
that  hydrogen  adatom  recombination  from  the  surface  of  this  covalent 
solid  occurs  via  the  irreversible  excitation  of  an  adatom  to  a  deloca¬ 
lized,  two-dimensional  band  state,  followed  by  reaction  of  this  deloca¬ 
lized  adatom  with  a  second  localized  hydrogen  adatom  to  form  a  hydrogen 
molecule  that  desorbs. 

Most  investigations  of  the  mechanisms  of  heterogeneous  reactions 
involve  studies  on  well-defined  surfaces  of  metal  single  crystals.  Only 
recently,  have  reactions  been  studied  in  comparable  detail  on  surfaces 
of  covalent  solids.  Silicon  surfaces  are  currently  being  investigated 
intensively  in  this  regard  [1-3].  The  Si(100)  surface  is  particularly 
appropriate  for  such  studies  because  the  stable,  reconstructed  (2x1) 
structure  is  rather  well  characterized.  In  the  absence  of  reconstruc¬ 
tion,  the  silicon  atoms  would  form  a  (lxl)  surface  structure  charac¬ 
teristic  of  the  bulk  lattice,  and  each  surface  silicon  atom  would  have 
two  "dangling"  bonds.  On  the  reconstructed  surface,  silicon  atoms  in 
adjacent  rows  "pair  up"  to  form  a  (2x1)  array  of  dimers,  resulting  in 
only  one  dangling  bond  on  each  surface  silicon  atom  [4].  In  the  mono¬ 
hydride  phase  of  hydrogen  on  Si (100),  there  Is  one  hydrogen  atom  bonded 
to  each  surface  silicon  atom  resulting  in  the  saturation  of  each 
dangling  bond  at  the  surface  and  the  preservation  of  the  (2x1) 


reconstruction  [5].  At  higher  coverages,  some  of  the  Si-Si  dimer  bonds 
are  cleaved  to  form  a  dihydride  phase  in  which,  at  saturation  coverage, 
one-third  of  the  surface  silicon  atoms  are  bonded  to  two  hydrogen  atoms 
[6].  In  the  studies  reported  here,  the  surface  coverage  was  suf¬ 
ficiently  low  that  the  recombination  kinetics  of  only  the  monohydride 
phase  were  studied. 

The  surface  coverage  of  hydrogen  during  isothermal  desorption  was 
monitored  as  a  function  of  time  via  laser-induced  thermal  desorption 
(LITD)  measurements  [7,8].  Briefly,  the  308  nm  output  of  a  XeCl  excimer 
laser  (~  2.0  raj  and  0^  »  75°  with  respect  to  the  surface  normal)  serves 
to  heat  rapidly  a  2.5  mm^  elliptical  spot  on  the  20  x  20  mm  Si(100)  sur¬ 
face  and  desorb  the  hydrogen  from  that  spot.  The  desorbed  molecular 
hydrogen  Is  detected  with  an  apertured  (17.5  mm  diameter)  and  differen¬ 
tially  pumped  quadrupole  mass  spectrometer.  Care  was  taken  to  insure 
that  the  mass  spectrometrlc  signal  from  a  laser  shot  scaled  linearly 
with  the  surface  coverage.  By  holding  the  surface  temperature  constant 
and  using  the  laser  to  desorb  the  hydrogen  from  a  series  of  spots,  the 
change  in  the  hydrogen  coverage  with  time  could  be  monitored.  Mo  damage 
of  the  surface  attributable  to  the  2.0  mJ  irradiation  has  been  observed 
in  over  six  months  of  repeated  studies.  Additional  experimental  details 
will  be  reported  elsewhere  [9]. 

The  apparatus  used  in  these  experiments  is  a  standard  stainless 
steel  bell  jar  30  cm  in  diameter  with  a  base  pressure  below  2x10“^® 

Torr.  Surface  order  and  cleanliness  were  confirmed  by  low-energy 
electron  diffraction  and  Auger  electron  spectroscopy,  respectively.  A 
tungsten  filament,  3  cm  from  the  silicon  surface  and  heated  to  1640  K, 
was  used  to  dissociate  molecular  hydrogen  (P^.  *  10“^  Torr)  during 


adsorption.  Thermal  desorption  mass  spectrometry  was  used  to  calibrate 
the  coverage  as  a  function  of  exposure.  The  desorption  peak  maxima  were 
at  783  K  for  the  monohydride  and  675  K  for  the  dihydride  with  a  heating 
rate  of  2.3  K/s.  In  the  isothermal  desorption  measurements  reported 
here,  the  surface  at  130  K  was  exposed  to  either  atomic  hydrogen  or 
deuterium  until  the  desired  initial  coverage  was  obtained.  The  surface 
was  then  heated  (18.3  K/s )  to  the  desired  temperature  (constant  to  ±0.2 
K),  and  the  LITD  measurements  were  made  as  a  function  of  time.  Initial 
coverages  of  1.0,  0.34  and  0.06  of  the  saturated  monohydride  were 
studied.  The  LITD  signals  were  normalized  to  that  of  the  saturated 
monohydride  coverage. 

The  decrease  in  the  fractional  coverage  of  chemisorbed  deuterium  is 
shown  as  a  function  of  time,  parametric  in  temperature,  in  Fig.  1,  for 
an  initial  coverage  of  0.34  of  the  saturated  raonohydride.  Similar  data 
were  recorded  for  both  hydrogen  and  deuterium  at  nine  different  surface 
temperatures  between  685  and  740  K.  The  data  were  quite  reproducible 
from  experiment  to  experiment  at  the  same  surface  temperature,  and  for 
different  initial  surface  coverages. 

If  the  desorption  reaction  were  second  order,  as  commonly  assumed, 


then  the  time  derivative  of  the  fractional  surface  coverage  of  hydrogen 
would  be  given  by  d%/dt  *  an<*  a  plQt  as  a  function  of 

time  would  be  linear.  Here,  kj  is  the  rate  coefficient  of  desorption  of 


hydrogen,  which  has  the  form  k<j  ■  Vd^exp^Ej/kgT).  On  the  other  hand, 

d0H 

if  the  desorption  reaction  were  first  order,  then  — rr  ■  ~  anc*  a 


plot  of  in  9^(0  as  a  function  of  time  would  be  linear.  The  data  were 
analyzed  in  this  fashion  and  a  typical  result,  D2  desorption  at  735  K, 
is  shown  in  Fig.  2.  Data  for  each  surface  temperature  and  initial 


coverage,  and  for  both  D2  and  H2  showed  similarly  good  agreement  with 
first-order  kinetics  and  similarly  bad  agreement  with  second-order  kine¬ 
tics. 

This  demonstrates  unequivocally  that  hydrogen  desorption  from  the 
monohydride  phase  on  Si(100)  Is  first  order.  All  traditional  notions 
concerning  adatom  "hopping"  on  a  surface,  followed  by  recombination  and 
desorption,  are  not  appropriate,  because  this  gives  rise  to  second-order 
desorption  kinetics.  There  is  a  completely  different  mechanism  for 
desorption  from  this  covalent  surface  compared  to  desorption  from 
metals.  Recent  isothermal  measurements  of  hydrogen  desorption  from 
Si(lll),  similar  to  those  shown  in  Fig.  1,  have  been  interpreted  in 
terms  of  second-order  desorption  kinetics  [8].  In  those  experiments, 
however,  the  fractional  surface  coverage  of  hydrogen  was  always  greater 
than  0.2,  which  corresponds  to  times  shorter  than  300  s  in  Fig.  2.  An 
analysis  for  this  limited  range  of  coverage  is  incapable  of 
distinguishing  between  first-  and  second-order  kinetics,  and  the  desorp¬ 
tion  reaction  is  likely  to  be  first  order  from  Si(lll)  as  well.  The 
fractional  coverage  range  displayed  in  Fig.  2  is  0.015  <  <  1,  and 

this  wide  range  is  necessary  and  sufficient  for  a  meaningful  analysis  of 
the  data.  The  fact  that  the  desorption  reaction  is  first  order  on  both 
surfaces  of  silicon  is  suggested  by  the  peak  shape  and  the  lack  of  a 
downshift  in  desorption  peak  temperature  with  increasing  coverage  in 
previous  thermal  desorption  spectra  [8,10-12],  but  most  previous  analy¬ 
ses  have  incorrectly  assumed  a  3econd-order  reaction. 

The  slope  of  the  straight  line  in  Fig.  2  is  the  rate  coefficient 
for  desorption  of  D2  at  735  K.  Similar  rate  coefficients  were  obtained 
at  each  temperature  for  both  H2  and  D2  desorption,  and  a  plot  of  in  kj 


as  a  function  of  reciprocal  temperature  is  shown  in  Fig.  3.  Based  on 
Fig.  3,  the  apparent  activation  energies  of  H2  and  D2  are  each  47  ±  1 
kcal/mol,  and  the  preexponential  factors  of  the  desorption  rate  coef¬ 
ficients  of  H2  and  D2  are  7.9  x  10^  s-*  and  5.6  x  10^  s”*,  respec¬ 
tively.  The  ratio  of  the  rates  shows  a  normal  kinetic  isotopic  effect, 
i.e.  k^OO/k^D)  ~  1.4.  It  is  also  important  to  note  that  the  fact 
that  all  constructions  such  as  that  shown  in  Fig.  2  are  linear  over  the 
entire  range  of  fractional  surface  coverages  (0.015  <  ®h(d)  <  D  implies 
that  the  rate  coefficients  of  desorption  of  H2  and  D2  are  not  a  func¬ 
tion  of  surface  coverage  [13]. 

It  is  necessary  to  explain  the  desorption  reaction  in  terms  of  a 
mechanism  that  is  consistent  with  the  observation  of  first-order  kine¬ 
tics  and  an  apparent  activation  energy  of  approximately  47  kcal/mol. 

The  challenge  associated  with  the  first  issue  is  self-evident.  The 
challenge  associated  with  the  second  issue  concerns  the  fact  that  in  the 
usual  analysis  the  activation  energy  must  be  larger  than  twice  the  Si-H 
bond  energy  [D(Si-H)]  minus  the  H-H  bond  energy.  The  measured  value  of 
Ej  then  requires  D(Si-H)  to  be  less  than  76  kcal/mol.  However,  values 
of  D(Si-H)  in  various  silane  molecules  are  all  much  greater:  a  value  of 
approximately  90  kcal/mol  for  this  bond  energy  would  be  expected  [14], 
The  following  mechanistic  model  resolves  both  of  these  difficulties  in  a 
straightforward  and  physically  appealing  way.  Hydrogen  adatoms  that  are 
chemisorbed  in  their  ground  state  on  Si (100)  are  thermally  excited  lrre- 
vers ibly  into  a  two-dimensional,  delocalized  band  state.  This  excited, 
delocalized  adatom  (H*)  then  reacts  with  a  localized  hydrogen  adatom  to 
produce  H2t  which  desorbs,  i.e. 


The  concentration  of  H*  is  always  small  if  the  first  step  is  rate 
limiting,  and  this  implies  immediately  that  &jj*  3  kj/k2,  and  that  the 
rate  of  hydrogen  desorption  is  given  by  Rj  ■  kjO^.  The  reverse  of  Eq. 
(1),  quenching  of  H*  to  produce  a  localized  Si-H  bond,  is  slow  due  to 
momentum  conservation  constraints  analogous  to  those  that  are  effective 
in  "selective  adsorption"  scattering  [16].  This  mechanism  explains  the 
fact  that  the  desorption  reaction  is  first  order,  and  the  meaning  of  the 
low  activation  energy  of  the  desorption  reaction.  The  activation  energy 
of  47  kcal/mol  for  the  irreversible  excitation  to  H*  reflects  the 
extreme  corrugation  of  the  H-Si(lOO)  potential  in  the  plane  of  the  sur¬ 
face,  cf.,  the  inset  to  Fig.  4.  This  should  be  expected  for  surfaces  of 
covalent  solids  with  directed  dangling  bonds.  Note  that  if  D(Si-H)  ~  90 
kcal/mol,  this  implies  that  the  H*  is  bound  (perpendicular)  to  the  sur¬ 
face  by  approximately  43  kcal/mol.  The  delocalized  H*  state  is  located 
further  from  the  plane  of  the  silicon  surface  than  the  localized  Si-H 
species,  and  therefore  the  effective  lateral  corrugation  that  H* 
experiences  is  weak. 

The  fact  that  the  proposed  mechanism  is  plausible  and  consistent 
with  all  experimental  observations  is  insufficient,  however,  to  verify 
its  validity.  We  have  demonstrated  the  validity  of  this  mechanism  une¬ 
quivocally  in  the  following  way.  Atomic  hydrogen  was  exposed  to  a  sur¬ 
face  containing  various  subsaturation  coverages  of  deuterium  in  the 
monohydride  phase.  The  deuterium  was  adsorbed  at  120  K,  whereas  the 


hydrogen  was  exposed  to  the  surface  at  680  K  for  100  s.  In  all  cases, 
the  deuterium  was  observed  to  react  with  the  hydrogen  and  desorb,  even 
though  680  K  is  far  too  low  for  any  significant  desorption  to  occur  from 
the  ground  state  surface  adatoms  during  this  period  of  time.  For 
example,  if  a  surface,  on  which  2  0.5  in  the  monohydride  phase  is 
initially  present,  is  exposed  to  atomic  hydrogen  at  680  K  (an  exposure 
which  would  almost  saturate  the  monohydride  phase  if  the  surface  tem¬ 
perature  were  120  K),  then  45Z  of  the  deuterium  originally  present  on 
the  surface  is  removed,  and  the  total  final  fractional  coverage  of 
hydrogen  and  deuterium  is  less  than  0.7.  In  the  absence  of  adsorbing 
hydrogen,  the  deuterium  coverage  would  have  decreased  by  approximately 
1%  at  this  temperature  for  the  same  period  of  time.  The  explanation  for 
this  observation  is  simple,  and  it  establishes  the  validity  of  the  pro¬ 
posed  mechanism  of  desorption. 

If  the  H*  species  exists  in  the  mechanism  for  desorption,  it  will 
also  be  a  precursor  in  the  adsorption  of  hydrogen.  However,  at  elevated 
surface  temperatures,  the  reaction  H*  +  D  +  HD  is  observed  (vide  supra). 
Assuming  that  D(Si-H)  ■  90  kcal/mol  requires  the  activation  energy  for 
this  reaction  to  be  greater  than  28  kcal/mol  and  less  than  47  kcal/mol, 
which  implies  that  the  activation  energy  for  dissociative  adsorption  of 
H2  is  less  than  19  kcal/mol.  This  picture  is  summarized  in  the  poten¬ 
tial  energy  diagram  of  Fig.  4.  The  qualitative  importance  of  this 
diagram  is  not  affected  by  adopting  the  reasonable,  90  kcal/mole,  value 
for  D(Si-H).  Relative  to  an  H2  molecule  in  the  gas  phase  as  the  energy 
zero,  two  localized  hydrogen  adatoms  on  the  silicon  surface  lie  at  -75 
kcal/mol,  one  delocalized  adatom  and  one  localized  adatom  lie  at  -28 
kcal/mol,  and  two  delocalized  adatoms  lie  at  +  19  kcal/mol.  The  dotted 


cur^a  shows  qualitatively  the  barrier  which  H2(g)  must  overcome  in  order 
to  chemisorb  dissociatively. 

The  "reactive  adsorption"  measurement  also  demonstrates  that  the 
reaction  between  two  H*  species  is  not  required.  The  H*  +  H*  reaction 
would  be  unlikely  based  solely  on  arguments  concerning  the  expected  con¬ 
centration  of  the  H*  species,  but  would  lead  to  a  rate  law,  R<j  * 
ki9jj/2,  which  is  indistinguishable  from  that  of  H*  +  H(a).  Hence,  it 
was  necessary  to  preclude  its  occurrence  in  order  to  establish  unam¬ 
biguously  our  proposed  mechanism. 

To  summarize,  a  new  mechanism  for  the  desorption  of  molecular 
hydrogen  from  the  monohydride  phase  on  Si (100)  has  been  proposed  and 
verified.  The  observed  first-order  desorption  kinetics  are  due  to  the 
irreversible  excitation  of  a  hydrogen  adatom  into  a  delocalized,  two- 
dimensional  band  state  on  the  surface  with  the  desorption  reaction 
occurring  between  this  excited  hydrogen  adatom  and  a  second,  localized 
hydrogen  adatom.  Both  the  kinetics  and  the  mechanism  of  hydrogen 
desorption  from  the  monohydride  phase  on  Si(100)  are  now  understood.  A 
similar  mechanism  would  be  expected  for  hydrogen  desorption  from  other 
surfaces  of  silicon,  and  from  surfaces  of  other  covalent  solids. 
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Figure  Captions 

Variation  In  fractional  surface  coverage  with  time  for  the 
Isothermal  desorption  of  D2  from  Si (100). 

A  comparison  of  first-order  and  second-order  kinetics  of  D2 
desorption  from  Sl(100)  at  735  K,  demonstrating  that  the 
desorption  reaction  Is  first  order. 

Arrhenius  plots  for  the  desorption  of  H2  and  D2  from  SI (100). 
The  implied  activation  energies  and  preexponential  factors  of 
the  desorption  rate  coefficients  are  47  ±  1  kcal/mol  and  7.9 
x  10*1  s”1  for  H2,  and  47  ±  1  kcal/mol  and  5.6  x  10**  s”*  for 
D2« 

Potential  energy  diagram  for  the  interaction  of  hydrogen  with 


the  Si(100)  surface.  The  corrugation  in  the  potential 
parallel  to  the  surface  is  shown  in  the  inset. 
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